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Abstract 

The six-quark instanton induced 't Hooft interaction, which breaks the 
unwanted ?7a(1) symmetry of QCD, is a source of perturbative corrections 
to the leading order result formed by the four-quark forces with the Ul{3) x 
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to bosonize the model by the functional integral method. We concentrate 
our efforts on finding ways to integrate out the auxiliary bosonic variables. 
The functional integral over these variables cannot be evaluated exactly. We 
show that the modified stationary phase approach leads to a resummation 
within the perturbative series and calculate the integral in the "two-loop" 
approximation. The result is a correction to the effective mesonic Lagrangian 
which may be important for the low-energy spectrum and dynamics of the 
scalar and pseudoscalar nonets. 
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1. Introduction 



In the absence of a quantitative framework within QCD to deal with its 
large distance dynamics, the physics of hadrons is usually approached either 
through an effective field theory or through phenomenological parametriza- 
tions based on some simple ansatz with solid symmetry grounds. In the first 
case the theory is written in terms of mesonic degrees of freedom Qll^l- In 
the second case the existence of underlying multi-quark interactions can be 
assumed. In this way one takes into account the quark structure of mesons 
explicitly, which may yield very important information about the structure 
of the QCD vacuum. Although these interactions are not renormalizable 
(as compared with CHPT, which is renormalizable in the sense of effective 
field theories), it is often possible to derive useful results in such models by 
introducing an ultraviolet cutoff. 

Not much is known about the origin of multi-quark vertices. The semi- 
classical theory based on the QCD instanton vacuum [3| provides evidences 
in favor of these interactions: two or more quarks can scatter off the same 
instanton (or anti-instanton), certain correlations between quarks originate 
from averaging over their positions and orientations in color space, the result 
being an effective quark Lagrangian. Assuming a dominant role for quark 
zero modes in this scattering process one obtains 2A^j-quark interactions [Nf 
is the number of quark flavors), which are known as 't Hooft interactions 
Actually, an infinite number of multi-quark interactions, starting from the 
four-quark ones, has been found in the instanton-gas model considered be- 
yond the zero mode approximation, all of them being of the same importance 
0. 

On the other hand, accurate lattice measurements for the realistic QCD 
vacuum show a hierarchy between the gluon field correlators with a domi- 
nance of the lowest one [6^. A similar hierarchy of the multi-quark interac- 
tions can be triggered after averaging over gluon fields. If this is true, there 
is an apparent contradiction with the instanton-gas model. This point has 
been stressed in 

The hierarchy problem of multi-quark interactions can be addressed on 
pure phenomenological grounds. For that we suggest to simplify the task 
considering only the four- and six-quark interactions. Once one knows the 
Lagrangian the obvious question arises: does the system possess a stable 
vacuum state and does this state correspond to our phenomenological expec- 
tations? If hierarchy takes place this question is pertinent for the leading 
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four-quark interaction, because in this case the effective quark Lagrangian 
can be studied step by step in the hierarchy with the assumption that four- 
quark vertices are the most important ones. In the opposite case we must 
study the system as a whole to answer the question. For the best known and 
simplest example, to which we dedicate most of our attention here, the solu- 
tion may be found analytically. As a result one can obtain definite answers 
to the above questions with a convincing indication in favor of a hierarchy 
for the considered example. 

Our choice of model is not accidental. The importance of the four-fermion 
interactions has been recognized for many years, starting in the early sixties, 
when Nambu and Jona-Lasinio (NJL) [S] used it for studying dynamical 
breaking of chiral symmetry. Later on a modified form of this interaction, 
written in terms of quark fields, has been used to derive the QCD effective 
action at long distances P ITIH [TT|. 

The six-quark vertices contain additional information about the vacuum 
|12j . They break explicitly the Ua{^) symmetry and are the only source of 
OZI- violating effects ^3|- Taken together with the NJL interactions, the 't 
Hooft Lagrangian gives a good description of the pseudoscalar nonet, espe- 
cially the rj and rj' masses and mixing ^1 ^j. In this form the model has 
been widely explored at the mean- field level [T^ . 

Let us discuss the Lagrangian which we shall use in our analysis. On 
lines suggested by multicolor chromodynamics it can be argued [T7] that 
the Ua{^) anomaly vanishes in the large Nc limit, so that mesons come de- 
generate in mass nonets. Hence the leading order (in Nc counting) mesonic 
Lagrangian and the corresponding underlying quark Lagrangian must inherit 
the Ul{S)xUr{3) chiral symmetry of massless QCD. In accordance with these 
expectations the Ul{3) x f//j(3) symmetric NJL interactions. 



G 



NJL 



{(iKq) + {qh^Kq) , (i) 



can be used to specify the corresponding local part of the effective quark 
Lagrangian in channels with quantum numbers = 0+, 0^. The Gell-Mann 
matrices acting in flavor space, Aq, a = 0, 1, ... ,8, obey the basic property 

The 't Hooft determinantal interactions are described by the Lagrangian 

in 

£h = «:(det qPLq + det qPRq) (2) 
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where the matrices Pl,r = (1 =F 75)/2 are projectors and the determinant is 
over flavor indices. 

The couphng constant k is a dimensional parameter {[k] = GeV"^) with 
the large Nc asymptotics k ~ 1/Nc^. The coupling G, [G] = GeV~^, 
counts as G ^ 1/Nc and, therefore, the Lagrangian (P) dominates over £h 
at large N^. It differs from the counting G ~ > which one obtains in the 
instanton-gas vacuum [Hj. 

It is assumed here for simplicity that interactions between quarks can be 
taken in the long wavelength limit where they are effectively local. The 't 
Hooft-type ansatz (j21) is a frequently used approximation. Even in this essen- 
tially simplified form the determinantal interaction has all basic ingredients 
to describe the dynamical symmetry breaking of the hadronic vacuum and 
explicitly breaks the axial U\{1) symmetry (TH]. The effective mesonic La- 
grangian, corresponding to the non-local determinantal interaction, has been 
found in [19J. 

Anticipating our result, we would like to note that if the hierarchy of 
multi-quark interactions really occurs in nature, the perturbative treatment 
seems adequate. The NJL interaction alone has a stable vacuum state corre- 
sponding to spontaneously broken chiral symmetry. But, as we shall show, 
the effective quark theory based on the Lagrangian^ 

C = qii^df, - rh)q + ^njl + -Cr (3) 

has a fatal flaw: if £h is comparable with £njl, it has no stable ground state. 
This feature of the model is invisible in a perturbative approach in £h- The 
situation is exactly analogous to the problem of a harmonic oscillator per- 
turbed by an term. This system has no ground state, but perturbation 
theory around a local minimum does not know this. In eq.Q the current 
quark mass, rh, is a diagonal matrix with elements diag(m„, m^, m^), which 
explicitly breaks the global chiral SUii^'i) x SUr(?>) symmetry of the La- 
grangian. 

There is also a special problem related with the bosonization of multi- 
quark interactions. To bosonize the theory one introduces auxiliary bosonic 
variables to render fermionic vertices bilinear in the quark fields. This proce- 
dure requires twice more bosonic degrees of freedom than necessary jT^j. Re- 
dundant variables must be integrated out and this integration is problematic 
as soon as one goes beyond the lowest order stationary phase approximation 

^The other multi-quark terms have been neglected here. 



4 



Pm OTj : the lowest order result is simply the value of the integrand taken 
at one definite stationary point ^H]- In this paper we shall show for the 
first time how to extract systematically the higher order corrections which 
contribute to the effective mesonic Lagrangian, and what to do with infini- 
ties contained in these corrections. The model is considered to illustrate 
our calculations. If the coupling k is not too small, corrections can be much 
larger than one might expect from l/Nf. counting, and be important for the 
mesonic (0"*", 0~) mass spectra and their dynamics. The large value of the 
77 — 77' mass difference obtained already at leading order and totally deter- 
mined by the six-quark interactions gives some credit to large corrections. 
Additionally, one can expect some enhancement at next to leading order by 
virtue of divergent factors: the cutoff scale is not known beforehand, and can 
be relatively large. 

Our study represents a very simplified view on the matter and should be 
considered as a first rough estimate which can be improved if necessary. The 
most essential approximation made here is related with the local character of 
the multi-quark vertices in eq.Q. It leads to 5(0) -singularities beyond the 
mean-field framework (see also Nevertheless, if it turns out that the 

low-energy QCD vacuum contains the hierarchy of multi-quark interactions, 
our approach can be used as a basis for a more serious work in this direction. 

The paper is organized as follows. In Section 2 the functional integral 
representation, for a bosonized version of the NJL model with 't Hooft 
interactions is derived. To make clear the approximations which will be used 
for the functional integration, a similar non functional integral, /, is consid- 
ered in Section 3. The reader who just wants to get a general idea of what 
we intend to study, can find it here. This section contains several clarifying 
points which are important in the following. In subsection 3.1 we give the 
exact result for /. Its stationary phase asymptotics is obtained in subsection 
3.2. The semi-classical method yields the same result, as shown in subsec- 
tion 3.3, as does the method presented in subsection 3.4. The point of these 
numerous calculations is just to show that all methods considered lead to 
the same asymptotics for the integral / which is totally determined by the 
number of stationary points. The perturbative treatment is given in subsec- 
tion 3.5. Here we discuss why the perturbation theory result is essentially 
different from the asymptotics obtained in the previous subsections, even for 
very small value of the expansion parameter. We resum the perturbative 
series in subsection 3.6. Calculating the next to leading order correction we 
illustrate the main goal of our studies in the forthcoming sections. 
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In Section 4 (subsection 4.1) we show that the chiral symmetry group 
imposes strong constraints which are only compatible either with the pertur- 
bative approach, or the expansion in a parameter that multiplies the total 
Lagrangian density (the loop expansion). Otherwise, as the consistent sta- 
tionary phase treatment shows (subsection 4.2), the model is unstable. 

The first alternative is considered in Section 5. To evaluate the functional 
integral we consider it like the natural infinite-dimensional limit of ordi- 
nary finite-dimensional Gaussian integrals (subsection 5.1). The perturbative 
treatment essentially simplifies calculations. Nevertheless, the integration 
over auxiliary variables leads to a special problem with 5(0) -singularities. 
We discuss this aspect of the bosonization in subsection 5.2. Another way to 
obtain the result is shown in subsection 5.3. Here a more elaborate spectral 
representation method is used to justify our computations. This subsection 
contains the prescriptions for the regularization of the 5(0) infinities, and a 
discussion of their reliability. 

The second alternative (the loop expansion) is considered in Section 6. 
Here, in subsection 6.1, we obtain in closed form the two-loop contributions 
to the functional integral Z and give arguments to justify this result. In 
subsection 6.2 we end with short conclusions and suggest future applications 
of our result. 

The summary is given and outlook is surveyed in Section 7. 
The Appendix contains the definition and the main properties of Airy's 
functions. 



The many-fermion vertices of Lagrangian Q can be linearized by introducing 
the functional unity [1^] 



2. Bosonization 



1 





a 



X 



exp <^ i / d''x[aa(sa - q\aq) + <\>a{Va - qil^Kq)] 



(4) 
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in the vacuum-to- vacuum amplitude 

Z = J VqVq exp J d^^xCj . (5) 

We consider the theory of quark fields in four- dimensional Minkowski space. 
It is assumed that the quark fields have color [N^ = 3) and flavor {Nf = 3) 
indices which range over the set i = 1,2, 3. The auxiliary bosonic flelds, 
(Ta, and, (pa, (a = 0, 1, . . . , 8) become the composite scalar and pseudoscalar 
mesons and the auxiliary flelds, Sa, and, pa, must be integrated out. 

By means of the simple trick (^, it is easy to write down the amplitude 
(0) as 

Z = [ VqVq n Vsa H T^Pa R ^^^a n exp U j d^xC) (6) 

a=0 a=0 a=0 a=0 ^ ^ 



with 

C = qih^d^, - a - i'y5(p)q + Sa{(Ta - rha) + Pa4>a + -C'njl + C-'n , (7) 

-^^NJL = f [{Sa? + (Pa)'] , (8) 

'^'n = [det(s + ip) + det(s - ip)] = —Aabc-Sa (sbSc - ^PbPc) ■ (9) 
d4 32 

We assume here that cr = aa^^a, and so on for all auxiliary flelds a,(j),s,p. 

The totally symmetric constants Aabc are related to the flavor determinant, 

and equal to 

Aabc ijfc^mni ('^a)jm('^b) j'n (-^c) fci 

= -^dabc + y 3 (S^aO^feO^cO " ^aO^ftc - ^Q^ac - ^cQ^ab) ■ (10) 

We use the standard deflnitions for antisymmetric fabc and symmetric dabc 
structure constants of f/(3) flavor symmetry. One can flnd, for instance, the 
following useful relations 

feacAbfc + febcAfac + fefcAabc = 0, 

deacAbfc + debcAfac + defcAabc = V^f^eO^afe/, 

8/^8 g 

^Aabb = -2W- 6aO, X! AaceAbce = " 5ab- (H) 

6=0 V c,e=0 y 



At this stage it is easy to rewrite eq.®, by changing the order of integra- 
tions, in a form appropriate to accomphsh the bosonization, i.e., to calculate 
the integrals over quark fields and integrate out from Z the unphysical part 
associated with the auxiliary bosonic variables (sa, pa) 



^ = / Y[l^(^a1^(t>a1^q1^q exp (^i J d'^xCq{q,q,a,^ 

X J Y[VsaVpa exp (^t J d^x£,(a,0,s,p)) (12) 



where 



A = qiil'^d^ -a- i'j5(l))q, (13) 

The Fermi fields enter the action bilinearly, thus one can always integrate 
over them, since one deals with a Gaussian integral. One should also shift the 
scalar fields aa{x) <ya{x) +ma by demanding that the vacuum expectation 
values of the shifted fields vanish (0|(Ta(x)|0) = 0. In other words, all tadpole 
graphs in the end should sum to zero, giving us the gap equation to fix the 
constituent quark masses rria corresponding to the physical vacuum state. 
The functional integrals over Sa and pa 

+00 

Z[cr, 0; A] = A/" y YiVsaVpa exp (^i J d^x£r(o" + rn, 0, s,p)j (15) 

—00 ^ 

are the main subject of our study. We put here = nia — rha, and Af is 
chosen so that Z[0, 0; A] = 1. 

Let us join the auxiliary bosonic variables in one 18-component object 
Ra = {Ra, Ra) where we identify Ra = Sa and Ra = Pa] a, a run from to 8 
independently. It is clear then, that R\ = + p^. Analogously, we will use 
IIa = {cTa, (pa) for external fields and Aa = (A^, 0). 

Next, consider the sum ^abcRaRbRc- If we require 

3 3 

^abc = —Aabc, '^abc = -J^^abc, $afec = 0, = 0, (16) 



we find after some algebra 



^ABcRaRbRc = ^^H (17) 
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with the following important property to be fulfilled 

<^abcSbc = 0. (18) 

Now it is easy to see that the functional integral (jl5j) can be written in a 
compact way 

+00 

Z[n, A]=M J 1[VRa exp (^i J d'^xCriU, A; R)^ , (19) 

-00 ^ 

where 

Cr = Ra{IIa + Aa) + jR\ + ^ ^abcRaRbRc. (20) 

We have arrived at a functional integral with a cubic polynomial in the 
exponent. 



3. Digression to a one dimensional case 

To get a rough idea of how to evaluate the integral ()19p we start with its 
one-dimensional analog 



+00 

/ = 



j dx e^-P(^) (21) 



where 

P{x) = ax + + Aa;3 ^ (22) 

and 0", a, h are constants. This integral plays the same role as our desired 
functional one, but is well defined as an improper Riemann integral of the 
real variable x. 



3.1 The exact result 



The integral ^21\ can be evaluated exactly. To show this let us express the 
polynomial P{x) in the form 

P"'(x ) 

Fix) = P(xo) + P'(xo)(x - xo) + -^r^i^ - ^^f , (23) 
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where Xq is chosen to satisfy the following equation 



P'\x) = bx + a = 

The coefficients are 

P{xo) -- 

Hence the integral (j2H) is given by 



P"'(x) 



+ 00 |- / 

/ = e'^^'"'^ J dx exp I Ip'{xo)x + 



3! 



-X 



We can rewrite eq. (j26j) in terms of the new variable t 



1/3 



X 



\P"'ixo)\, 



thus we arrive at 



I = e*^(^o) I ^ 



1/3 



dt exp ( M "2" '^^ 



where 6 is defined to be 



5 = sgn(6) P'(xo) 



1/3 



(24) 



(25) 



(26) 



(27) 



(28) 



(29) 



This integral is well known. The result of integration can be represented in 
terms of the Airy function (see Appendix for details) 



/ 2 \ 

/ = 27re*^("°M— j Ai{6). 



(30) 



If 6 is real, the function Ai{6) is also real, and the phase of / is equal to P{xo). 
In the following, after some generalizations made for the integral /, its phase 
will represent the effective action of a dynamical system. The expression for 
this phase is the main goal of our calculations. 
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3.2 The stationary phase result 



The result (jHUj) for / is exact. It can be approximated at large values of 5 
by its asymptotic series, still giving us an exact expression for the phase. To 
obtain asymptotics let us transform eq. ()28|) to a more convenient form by 
the replacement 

t = My. (31) 



One has 




1/3 



+ 00 

\5\ I dy exp {i\5\ 



3/2 



9{y) 



(32) 



where g{y) = y^/3 + sgn.{6)y. The integral (jH^ is already in the form appro- 
priate for the stationary phase method to be applied. Indeed, as S —oo 
the term \5\'^^'^g{y) gives a rapidly oscillating contribution to the integrand 
in eq. ()32j) which cancels out, except at the regions of critical points. The 
contribution from these regions can be evaluated on the basis of the station- 
ary phase method. In our case a ~ G ~ & ~ k ~ ^/N^, cr ~ 1, i.e., 
^3/2 ^ jys large Nc arguments^ can be used to justify the 6 ^ —oo 
limit. Both critical points {g'{y) = ^ yi^2 = ±1) belong to the interval of 
integration. Thus we have 



[\b\. 



1/3 



-t- exp 



\6\ J dylexp 




(5 



-oo). 



(33) 



The last term in both exponents can be factorized and then expanded in a 
power series of 



1/3 



+ 00 



\b\. 



\S\ / dyiexp 



+ exp 



_,|5|3/2^,2_- 



y 



}£;^(#rfl" (^---). (34) 



^AU large Nc arguments given in this section hint of course at the functional integral 
case which will be discussed later. 
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and integrated term by term. The corresponding integrals are evaluated 
exactly 



/ Ay exp {±^ay') y'^ = (±^)"^^j^^ e^^/^ (^) (a > 0) (35) 

— OO 

giving as a result the following asymptotics for the integral / 



|<5|l/4|6|l/3' 



cosgi5r-^)+o(i5r^/^) 



(5^ -oo). 



(36) 

Let us stress that both critical points contribute to the resulting asymptotic 
series reproducing the well known asymptotics of the Airy's function ()156|) . 
We have considered here the case when the series oscillates. In the opposite 
case, 6 oo, the asymptotics falls down exponentially (jl55|) . and does not 
have the appropriate form from the physical point of view. 

3.3 Semi-classical asymptotics 

One can calculate the integral ()2H) assuming that a large parameter is already 
present in the exponent. For instance, the reduced Planck constant, h, can 
be considered as such a parameter 

+ 00 

I = J dx exp Q^(2;)) ■ (37) 

— oo 

In this case one obtains the asymptotic expansion for I at small values of 
h 0. This approach is known as the semi-classical expansion. 
The real function P{x) has two critical points xj {j = 1,2) 

P'[x) = -x^ + ax + a = => Xi2 = —, • (38) 

2 ' b 

Both of them are real at > 2ba and, therefore, belong to the contour 
of integration. In the neibourhood of these points the polynomial P{x) is 
conveniently written in the form 

P"(r) P"'(r) 
P{x) = P{x,) + -^{x - x,r + ^-^(x - x, f , (39) 



12 



where 



\ h f-iV 



D = a^-2ha> 0, 



P"(xj) = (-ly+i/^i/^^ P"'{x) = b. 
Thus the integral under consideration is estimated as 



2 +?° r / p///^ 

I^Y^e'Pi-^^ J dxexp ^ ^ 
i=i L \ 



or 



DV4 



COS I -777^ - T + ^(^) 



0). 



(40) 



(41) 



(42) 
(43) 



V 362 4^ 

Noting that 

one can see that this result coincides with our previous estimate (jHBj) . 

It may be helpful to remark that if one would take the contribution of 
only one critical point in (j^T|) . the resulting asymptotics would be obviously 
different, and, what is important for us, the phase too. 



3.4 The Nc 00 asymptotics 

The representation ()39|) can be used as a first step to estimate the integral 
/ without any reference to the semi-classical expansion. Alternatively, our 
calculations can be based on the large N^. asymptotics. According to the 
formula (j39p . we have identically 



/ = e'^^^'^ / dx exp i 



I 3! 



This holds whether xi or X2 (see eq. fj38|) ) is chosen here. 
Next, replacing the variables 



X 



3P"{x,) 



P"'{x,) 



we arrive at 



3VD 



dt exp (iA/(t)) 



(44) 



(45) 



(46) 
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where we have used the following notations 

92^3/2 



A 



262 ' 



(47) 



This form of the integral is already appropriate to apply the stationary phase 
method at large A. Large N^. arguments can be used to justify the A ^ oo 
limit, for it is known that A ~ N^. 

Let us obtain the leading term in this asymptotics. The critical points 
are given by the equation 



f'{t)=0 ^ h = 0, t2 = {-iy sgnib)-. (48) 
Therefore, the integral / has the following asymptotical estimate at A ^ oo 



r— +00 

I ^ e^^(^^)^ / dt {exp [^A ((-1)^+42 + sgn(6) t^) 

\u\ J 



+ exp 



-ir^'^ + (-ir^^ + sgn(6)t=')]} 



(49) 



The first term of the integrand is the contribution of the critical point ii, the 
second one comes due to ^2- Noting, that 



one has 



P(a;,) = P(a;o) + (-l)^-^, 



(50) 



j^e^P(.o)^ \dt {exp \\ ((-1)^1 + {-irn-) 



+ exp 



-iA ((-1)''^ + (-1)'+'^')] } e^^^S^^*-)*' 



(51) 



It is now obvious that the integrand does not depend on j, thus 



iP(xo) 



I ^ e 



+ exp 



3/D 



-i\ 



27 



dt \ exp 



-t' 



ixl-- e 

. V27 



i\ Sgn(6) 



(52) 
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One can see finally that the obtained asymptotical series 



\b\ VA 

coincides with our previous result 



e-''"''4Fs^cos{----^]+... (53) 




3. 5 The perturbation theory approach 

Our asymptotic estimate requires a further explanation in what concerns the 
definition of the improper Riemann integral (j2ip . It is understood as the 
limit 

L 

/=lim/(L), I{L) = f dxe'^^""^ (54) 

where we assumed that both stationary points xi^2 (see eq. fl38p ) belong to 
the interval —L < x\^2 < L; the asymptotics of I{L) was formed accordingly 
by two independent contributions: I{L) ~ Ixi{L) + I^c^i^L). Integrating in 
Lj..{L), we took into account that only a neighbourhood of a stationary point 
Xi is important, and extended the integration limits to L ^ oo. 

Let us suppose now that the coupling h is small, as compared with a, and 
one can consider the limit 6 — >^ 0. It is easy to see that 

xi^-- + 0(6), X2~-^ + - + C(6) (55) 

at small b. The first solution is regular at 6 — > 0, but the second one shows 
the singular behavior 

hm IX2I = +00 . (56) 

This behavior refiects the increasing importance of the stationary point xi 
in the physical problem for which the cubic term ~ 6 is considered as a 
perturbation. The second stationary point, X2, varying as 1/b, can finally 
leave the interval {—L, L), changing as a result the asymptotics of the integral 
J(L) to 

/(L)~4,(L), X2^{-L,L). (57) 

The intuitive argumentation given above must be clarified. One should 
not think that the real value of L is very important for the matter. Actu- 
ally, one excludes the singular critical point from the phase P{x) by directly 
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expanding the phase in the neighbourhood of a regular stationary point, as 
we aheady did it for /^^. What is really important here is to conclude that 
only the regular critical point determines the perturbative regime of the sys- 
tem. According to this attitude, we use the term "perturbative regime" in 
a wide sense: the standard perturbative expansion in powers of b, which we 
will obtain below, is merely another way of looking at the one critical point 
asymptotics of I. 

Let us separate the unperturbed part Pq from the perturbation Ppt in the 
integral / 



P{x)^Po(x)+Ppt{x), Po{x) = ax + -x\ Ppt{x) = -x\ (58) 
Expanding the integrand in powers of the coupling b, we obtain 



n=0 



y-p 



n=0 



pt 



exp 



ib 
6" 



. d ' 



+00 



+ 00 



(59) 



The new feature is the fact that the integrand of /pt has only one stationary 
point, Xst = —a/a, as opposed to /. The singular critical point is gone. 
The evaluation of the Gaussian integral is straightforward 



+ 00 




27r / . / TT a 

exp t 



dxe^-P°(^) = 

) 

This reduces our integral Ipt to the form 



2^ 

2a, 



'pt 




a 



ib / . d 
da , 



6 



—I 



da , 



+ ... 



exp 



72 

'2a, 



(a > 0). 



—I 



da 



(60) 



(61) 
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Up to the terms of order we have now 

r3 \ / 



+¥ 



, , a ta" \ I 5i 5cr 5zcr 



35icr 35a^ 7icr^ 



+ 



a 



24a3 8a^ 2Aa^ 72a^ 



9 \ 



(4805 ' 36a6 24a7 3Qa^ ' 1296a9 



1 + 



(bai + 6^a2 + &^«3 + • • •) 
1 + b{/3i + iai) + 6^ (^i{a2 + aiA) + /^a - y j 



Noting, that 

(1 + bPi + b^P2 + b^p3 + ...) exp 



(62) 



+ b^ 



i {a3 + ai(32 + a2/5i - -77- ) + - «i«2 



6 



y-i 



1 2 
a. 



+ ... 



one can represent this result in the more convenient form (up to the same 
order of accuracy) 



/pt = e^-/^ 




27r / ^ a ,2 . ISa^ 



3«-Poff(o-) 



2a2 



16a6 



(63) 



where 



b — - + b^ 

2a 6a3 



0" 



5cr a 



24a3 8a 



^ +^^hr:^-^ +W. (64) 



The perturbative result Pefr(o") is an approximation, which can be sys- 
tematically improved. In its truncated form the obtained series differs from 
the expansion which one can obtain from I^^ (see next subsection for de- 
tails). It is clear, however, that if one sums up the series, the perturbative 
and asymptotic estimates will coincide, i.e., /pt = Ix^. 



3.6 Resumming the perturbative series 

Let us calculate the contribution to the integral / which comes from the 
critical point = 1 in eq. ()34j) . It is not difficult to find out that exactly this 
point represents the regular solution at 5 ^ 0. For this part of the integral 
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/ we will use the symbol /„j , as we did before for /a- - , 



T e 



iP{xo) 



1/3 



\5\ I dy exp|^|5|3/2 (y' 



(65) 



n=0 



It is clear that in the last sum only the even powers of n contribute. For 
the first two terms (n = 0, 2) one obtains 



1+ .JL,. +--] ■ (66) 



48|5|3/2 



Since 



1^1 



3/2 



262 



1 - 



2ba 



3/2 



a2 ; 



we can rewrite this result in the compact form 



(67) 



ly, ~ e 



in /A 




A{a) e<« ( 1 + 



hi 



48|5|3/2 



(68) 



where 



A{a) 



2ha 



-1/4 



1 + 



2^^~8^ 



+ 



156V3 
16a6 



(69) 







|-( 




362 




6a3 


6V4 


2a ~ 


6a3 


8a5 



26a y/^' 
a2 j 



63^5 



+ 0(6^ 



(70) 



One can explicitly see that the series in powers of 6 obtained from the integral 
ly^ perfectly coincides'^ with the perturbative expansions (jFK?jl and (jMj) . 



The series fl68|l corresponds to a partial resummation of the perturbative 
series. In order to see what is resummed, let us introduce a parameter t 



^Actually, we checked the terms up to and including It' order. 
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according to the substitution: a,b,a ta,tb,ta in eqs. ffH3|) and We 
find for Jpt /pt(t) 

^ eW. ^ (l H- . . 5^ . i|! ^ . . .) e-(-., (71, 



wliere 

^ , , fa^ ha^ b^a' b'a'\ 1 / 5b' 56V\ ,^,,4, 
^eff(t, ^) = - k7 + 7r::T + ir:r + ir:F +77 7773 + 177^ +(^(&')- (72) 



2a 6a3 Sa^ Sa^ / t2 \^24a3 



6a 



It is clear now tliat tfie leading term, Pg^, is represented in Pes(t) by tfie sum 
of terms which are oc the first stationary phase correction is given by 

terms of order oc and so on. The field theoretical analog of this power 
series expansion in (l/t) is well known as a loop expansion [22]. 



4. Outlining the problem 

So far we have worked with a one-dimensional case. We now want to clarify 
what actually can be learned from this simple example for the functional 
integral (fTn|) considered. 

It is clear, for instance, that we cannot follow the approaches discussed 
in subsections 3.1 and 3.2. An analog of eq.()24|) for our functional integral is 
the equation 

^ C^'ab = G5ab + ^^ABcRc = (73) 

which cannot be fulfilled. 

On the contrary, the system of equations based on the first order deriva- 
tives 

-^ = GRa + Aa + Ua + - ^abcRbRc = (74) 

is self-consistent and can be solved Therefore, one can obtain the semi- 
classical asymptotics by analogy with the stationary phase method discussed 
in subsection 3.3. It corresponds to the case without hierarchy, since two 
critical points are considered. We will show, however, that such a system does 
not have a stable vacuum state, i.e., the hierarchy of multi-quark interactions 
is a very important physical requirement for the model. Let us discuss the 
matter in detail here. 
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4-1 Solving equation fl^ 

We need to recall shortly the solutions of eq. lff^ . Up to some order in the 
external mesonic fields, 11^, we may write them as a polynomial Ra = TI^a 

7^? = Hf + H%IiB + H%c^b1^c + H^scd'^bI^c^d + • • • (75) 

where i = 1,2, . . . denote different possible solutions. The coefficients H^^ 
depend on and on the coupling constants G, and the higher index 
coefficients H^^ are recurrently expressed in terms of the lower ones. For 
instance, we have 

{h^b) = - {G5ab + f^^^ABcH^^) , (76) 
Habc = 2^defH^daH^ebHfc > (^7) 

and so on. 

Putting this expansion in eq. ()74p one obtains a series of self-consistent 
equations to determine H^^ . The first one is 

GH^^ + ^A + ^ ^abcH^^h'^ = 0. (78) 

One can always find the trivial solution Ha = 0, corresponding to the 
unbroken vacuum = 0. There are also non-trivial ones for the scalar 
component, i.e., 

Hf = ih^\0). (79) 

The number of possible solutions, i, depends on the symmetry group. 
The coefficients /i^*^ are determined by the couplings G, k and the mean field 
value Afl. In accordance with the pattern of explicit symmetry breaking the 
mean field can have only three non-zero components at most with indices 
a = 0, 3, 8. If two of the three indices in Aabc are {0, 3, 8}, then the third one 
also belongs to this set. Thus A^ is the only object which determines the 
vector structure of the solution hf>, and therefore /i^*^ 7^ if a = 0, 3, 8. It 
means that in general we have a system of only three equations to determine 
/i« = /,«A„ = diag(/z«,/i?,/i«) 

GK + A, + ^hdh, = 
16 

Ghd + + ^KK = (80) 
16 

Ghs + As + —Khd = . 
16 
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This system is equivalent to a fifth order equation for a one-type variable 
which can be solved numerically. For two particular cases, when m„ = = 
rhs and m„ = rhd ^ rhg, eqs. ljHUj) can be solved analytically [21] • The simplest 
example: = rhd = t^s (or, equivalently, h^^ = h^^ = /i^*^) corresponds to 
SU{3) flavor symmetry. In this case eq.ffSUj) has two solutions 

If 46"^ > kA, they are real and will contribute to the stationary phase tra- 
jectory. 

4-2 The lowest order semi-classical asymptotics 

Since the system of equations (ff^ can be solved, we may replace variables 
Ra — > Ra = Ra — Tl^A the functional integral (fT^ to obtain the semi- 
classical asymptotics 

Z[n,A]~ArX: exp(^/dV4tO 
1=1 ^ ^ 

X j WVRa exp j d^xC%{n^^)RARB) 

-oo ^ 

X E ^ y d^xRARsRc) {h ^ 0) (82) 

where n is the number of solutions, 7^^"*, of eg. ((71)). has been defined in 
eq. lfTHj) . and 

= f (7^?)'+ ^ (n^ + Aa) = + OiU'). (83) 

Here we used eq. lff^ to eliminate the term proportional to k. Let us stress 
that depends on k implicitly: k, is contained in TZ^a ^ or more precisely in 
the coefficients H^^ which are functions of h^^\ This dependence is singular 
at K — > 0. One can see this, for instance, from eq. ()8H) where h^^ — > oo for 
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small K. This behavior reflects the fact that we are far from the perturbative 
regime, meaning that the interactions £njl and £h are equally weighted. 

The linear term in the a fleld is written explicitly. This part of the 
Lagrangian is responsible for the dynamical symmetry breaking in the multi- 
quark system and taken together with the corresponding part from the Gaus- 
sian integration over quark flelds in eq. lfT^ leads us to the gap equation. 

At leading order, = 0, we have the estimate for eq. (j82j) 



aa + ...] 



where y4(j) is real and proportional to 

A(,)~|det£:;^(7^«)|-l/^ (85) 
Therefore, if one considers the case with m„ = = rhs, Z is given by 

Z ~exp (^t J d^x^/i«a, + ...^ 

8G f \ 
I— J d^x{au + (rd + (rs) + ■■■] (86) 



exp 



Let us recall that the quark loop contribution to the gap equation is 
well known (see, for instance, ^22j|). Combining this known result with the 
estimate (|HHjl . one can obtain the corresponding effective potential U{m) as 
a function of the constituent quark mass m 




Uim) = m m A — m In 

+ ASln(l + ^)|, (87) 

where we consider the case m = for simplicity and A^ denotes the cutoff of 
quark loop integrals. 

We plot this function in flg.l for G/ k > 0. The parameters are flxed there 
as Aq = 1 GeV, G/ k = 2.34-10"'^ GeV'^. One sees that the system has at most 
a metastable vacuum state®. We must conclude that the model considered 
has a fatal flaw and can be used only in the framework of the perturbative 
approach, which does assume the hierarchy of multi-quark interactions. 

^In the opposite case, G/k < 0, the effective potential does not have extrema in the 
region m > 0. 
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Figure 1: The effective potential U{m) (see eq. (jH7|) ) for tlie parameter set 
given in the text as a function of the constituent quark mass m. All are given 
in dimensionless units. 

5. Perturbative expansion of Z 

We shall restrict ourselves in this section to the perturbative treatment of 
the functional integral |T9|l . The loop expansion will be considered in the 
next section. 

5.1 The perturbative series 

Let us devide the Lagrangian (pUj) in two parts. The free part, Cq, is given 
by 

£o(^a) = Y^A + ^A(nA + Aa). (88) 
The 't Hooft interaction is considered as a perturbation 

£j{Ra) = ^ ^abcRaRbRc- (89) 

Thus the perturbative representation for the functional integral (|T9|l can 
be written as 
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where 

= (91) 

Since the boson fields appear quadratically in eq. ()90p . they may be inte- 
grated out, yielding 

Z = Nexp(tJ Cj{X)^ exp (^-t J (92) 

where 11^ = Ha + A^. The overall factor = {—2'n-i/G)^N' is unimportant 
in the following. 

We want to calculate the effective action T^s, which by definition is the 
phase of Z 

Z = AiTlA) exp (^z J CesiTlA)^ , (93) 
and ^(11^) is a real function. Comparing (jU^ and one gets 

Teff = ^ In ^ + To - z In (1 + e'^^o [e' I^' - l) e*^") . (94) 
Here Tq represents the leading order result for Fcfr 

1 r^. 



ro = -^ynl (95) 

while the second logarithm in eq. ()94|l is a source of f/A(l) breaking corrections 
which arise as a series in powers of the functional derivatives operator 

ri = J Ci{Xa) = I ^ABC J XaXbXc . (96) 

To make this statement more explicit let us consider the expansion 

00 -m 

S = e-^ro (e^r, _ ^^ro = g i- (e-^^T^e^^") . (97) 



m=l 



Taking into account the symmetry properties of the coefficients ^abc and 
our previous result (fTHj). we find 

+ I UaXbXc - XaXbXc^ (98) 
24 



so that 

oo 

6 = '^"^n, (99) 

n=l 

represents the effective action ()94|) as a perturbative series in powers of k. 
For instance, up to and including the second order in k we have 

Teff = i In 4 + To - m5i - f ^2 - Idf) - ... (100) 



where 



N 



Si = iT77^ "^ABc I n^n^nc , (loi) 



3!G3 

The real factor A{Tl) is always chosen such as to cancel the imaginary 
part of the effective action. To the approximation considered we have, for 
instance. 



A{n)=NexpUo) (^)7n 



= ^(^i + ^(o)^ynl + ...j. (103) 

It contributes to the measure of the functional integral over cxa, 0a- 



5.2 The problem with infinities 

The terms with the 5(0) function require further explanation. The fact that 
auxiliary fields can lead to special problems with infinities is well-known 
pij . For instance, a relevant case can be found in pS]: the new variable A is 
introduced in the photon self-energy diagram in a similar way, i.e., through 
the integral over a 5-function, see eq.(8.18) of the book, which is equal to 1 
and finally, after changes of variables in the whole expression eq.(8.19), the 
logarithmically divergent integral over A requires a cutoff. 

A case involving a functional integral is given in [21] where the author 
considers as an example the non-linear cx-model and shows that the Gaussian 
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functional integral with a field dependent coefficient in the exponent (there 
are no derivatives of fields there) leads to a 5(0) factor in the corresponding 
effective action (see also P?]). 

In our case auxiliary variables have been introduced by eq.(@]) and chang- 
ing the order of integration in eq.()12|) we finally came to a similar problem 
with infinities. 

To see the origin of the encountered singularities in detail we shall use the 
language of Feynman diagrams. The graphs contributing to lowest order in 
K are shown in fig. 2. In these diagrams, a line segment (straight or curved) 




Figure 2: The lowest order graphs contributing to 6i. 
stands for a "propagator" 

Aab(x -y) = -{i/G)5ab5{x - y), (104) 

extracted from the last exponent in eq. ()92|) . A filled circle at one end of a line 
segment corresponds to the external field, i J d'^xIlA{x), and a vertex joining 
three line segments is used for in^ABC I d'^x. 

The first diagram represents the 6i term in eq. ljlOlj) . The contribution 
of the second tadpole diagram is equal to zero. Indeed, the vertex contains 
the group factor ^abc- The contraction of any two indices in this factor by 
6ab from the propagator (situation occuring for the tadpole graph) reduces 
it to zero, according to eq.lfTHj). We thus find that tadpole diagrams do not 
contribute due to the fiavor structure of the 't Hooft interaction. 

To next to leading order in k we have the four graphs shown in fig. 3. 
As we just learned, the third diagram does not contribute. The other ones 
correspond exactly to the three terms of eq. p02|) . The first tree diagram is 
finite. The second one-loop diagram has a divergent factor 5(0). The last 
two- loop diagram contributes as 5(0)^. These singularities were caused by 
the local structure of the last exponent in eq. (j92j) or, which is the same, by 
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Figure 3: The diagrams of the order contributing to 

the 6{x — y) term in the propagator ()104|1 . We beheve that if one would start 
from non-local NJL interactions, the singularities could be weaker or would 
even disappear. 

The factor 5(0) requires a regularization. This is an expected trouble in 
the NJL model which is nonrenormalizable and, as a consequence, the funda- 
mental interactions must be cut off. The cutoff is an effective, if crude, imple- 
mentation of the known short distance behavior of QCD within the model. 
In the next subsection we shall discuss this procedure in the framework of 
the spectral method. The value of the cutoff can be fixed by confronting the 
model with the experiment. In this interpretation the model is effectively 
finite, including the higher order corrections.^ 

The diagrams are a very convenient language to understand another fea- 
ture related to the multi-loop contributions: the phase factor corresponding 
to the diagram can be simply calculated. Indeed, it is easy to see that for 
any diagram the formula 

E = 2I -W (105) 

is fulfilled. Here E is the number of external fields, / stands for the number 
of internal lines, and V is the number of vertices. On the other hand, the 
number of loops, L, is given by 

L = I -V - E + l. (106) 

This is because the number of loops in a diagram is equal to the number of 5 
functions surviving after all integrations over coordinates are performed, ex- 

^Actually, if one works with the NJL model, one must choose among several known 
regularizations. Unfortunately, the dimensional regularization (DR) cannot be used. The 
gap equation in this case does not have solutions and as a result there is no dynamical 
chiral symmetry breaking. This is why we cannot simply take advantage of the well-known 
result (5(0) = (in DR) to avoid the problem. 
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ept for one over-all integration related to the effective action. Every internal 
line contributes one S function, but every vertex or external field carries an 
integration over the corresponding coordinate, and thus reduces the number 
of 6 functions by one. 

This result shows that the overall phase factor of a given graph i^-^-^ = 
is entirely determined by the number of loops. In particular, diagrams 
with an even number of loops contribute to the effective action (the argument 
of Z), while the diagrams with odd number of loops contribute to the measure 
(the modulus of Z), see eq. (j^ . 



5.3 The spectral representation method 

The problem with 5(0) singularities can be analysed with more accurate 
methods. In this subsection we shall use the spectral representation method 
to evaluate the integral (fT^ . 

Let us consider, for simplicity, a one-dimensional field theory approxi- 
mation for the integral (fT3j) . The number of dimensions does not matter 
here. The extension to the real case is straightforward. We also put our 
system in an interval of size L, i.e. —L/2 < x < L/2, assuming the limit 
L — s> oo in the end of calculations. The Fourier decomposition of the fields 
fa{x) = {sa,Pa,(^a,(pa} insidc the interval is 

faix) = E fn exp Un^) . (107) 

n=-oo ^ / 

This corresponds to the periodic boundary conditions /a (—L/2) = /a (L/2). 
The Fourier coefficients represent the field faix) inside the considered 
interval 

= 2J^^^ faix) exp (-^277— j , = (108) 



Then we have 



/ 



/■L/2 

dx Cix) - J ^^^dx Cix) = LC , (109) 



where C is devided in two parts: the free one Co, and the perturbation Cj 

(110) 



+ 00 p 

^0= E y(<^-n+KP-J+<^-„ + 0>-, 
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>C/ = ^A^bc E < - ^P'mPi) Sn+m+k,0. (Hi) 

n,m,k 

The functional integral can be understood as the product of integrals over 
the Fourier coefficients 

+00 +00 

f VSa{x)Vpa{x) ^ n / (112) 
—00 —00 

and in the perturbative approximation can be written as 

+00 

Z[a,0]^exp{^L£,(5;^,P„")}n / d<dp;; exp (^L£o) (113) 

—00 

where 

5" = -- ^ = -- ^ (114) 

It is convenient to use the normalized functional Z defined as 

(lis) 

z[o,o]- ^^^^^ 

The Gaussian functional integrals can be simply evaluated, for instance, 
one has 

-^M = n y exp |iL E y-^SnS-n + CTnS-nj j 

+00 ^ 

= J dso exp jiL (^—sl + aoSo^ | 
—00 

00 +?° r +00 ^ 

X [] y ds^ds;* exp <^ iL ^ (CszS;* + (^/s^ + a^si) \ 



iL 



+00 



= /[0]exp|-— J2 ^no-ny (116) 

Thus the functional integration yields 

{jj +00 ^ 
E (^x„a_„ + 0n0-„)|. (117) 
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We should now calculate partial derivatives to obtain the effective action 



Z ~ A[a, (p\ e 
As before (see subsection 3.2) we have 



Teff = i In A[a, 0] + Tq - z In (l + e-'^° (/^' - 1 



118) 



(119) 



where A[a, 0] is fixed by the requirement that the effective action T^ff is real. 
Here Fq represents the leading order action 



T OO 



(120) 



The logarithm in eq. ()119j) is a source of Ua{^) breaking corrections which 
arise as a series in powers of the partial derivatives 

rj = LCj{s:,p:) 

n,m,k 

One has the following result 



-ir, 



or J'^o — A \^ js 

■L/C — —-— /iabc On+m+k,0 
n,m,k 



G3 



my^k/ 



3 6 b ■ ^ b 

'q^ {^m^k ~ 4>m't>k)^n + ^ ^n4>mPk 



3 fi 

qaf qb qc O pb TDC 

'-'nWm'-'k '^■^m-^k 



Therefore, expanding 



(122) 



n=l 



one can determine 



Si 



-iL 
32G3 



^abc E ^n(0"icrfc - 30^0fc)5n+m+fc,O: 
n,m,k 



(123) 



(124) 
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Sf . 9L 
' 2 64(2^)5 



^abc^abc ^ ^m.+A;+m+A;,0 



(126) 



+ E + ^r.r^n) E 1 + ^ (e i) ^ (125) 

and so forth. 

The effective action can be also expanded in k 

where one assumes that ln74[cr, 0] must be also expanded, i.e., 
A[a, 0] = 1 + + t^l32 + 0{K^), 
In A[a, 0] = + I^^^ " f) + (^('^'), 
and, up to the considered order in k, we have 

A =0, /?2 = 7^ E (<^-n + C0-n) E 1 ' 

Therefore the perturbative action systematically obtains Ux{l) breaking cor- 
rections in k", which we collect in the corresponding part of the action, r„. 



(127) 



(128) 



hoo 



eff 



ETn, 



(129) 



where Fq is given by eq. (|12(jp . and 
Ti - 



A \ " „a(h a Q^fe j^c 



326-3 
"64(26) 



n.m.k 



5 AabcAg^l,^ ^ 5m_l_fc_l_^+fc,0 
m,fc,m,fc 



(130) 



Here an unessential constant (for the physical action) has been omitted in 

r2. 
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Using the well known relation 



exp ii——n\ = L ^ S{x — Ln) = LS{x) (131) 



where on the last step we used the fact that in the problem considered all 
x-dependent functions are integrated only in the interval —L/2 < x < L/2 
and, therefore, only the term with n = can contribute, one can obtain that 



E/n/tn = 7 / dxfa{x)Mx)exp{-t27rkx/L). (132) 
As a result, taking the infinite- volume limit L ^ oo, we have, for instance, 

^ +00 

— I dx {al{x) + <Pl{x)) . (133) 



2G 

Using the following relation 



E fnlQ ■ ■ ■ fZ '^ni+n.+...+n.O = jl dx r{x) ^x) . . . ^ (x) 
n2,...,ni ^ J-L/2 

(134) 



in (|13(J|) . and taking the limit L ^ oo one finds that 



Ti = ~^^^^a6c / dxaa{x) [ab{x)ac{x) - 3(f)bix)(f)cix)] 
AabcAabc I dx [Aab{x)ai{x)(j)c{x)4>c{x) 



64(2G)5 

+ {ab{x)ac{x) - (f)b{x)4>c{x)) {a-b{x)ac{x) - 05(x)0c(x)) ] . (135) 

The infinite- volume limit for the expansion coefficient /?2 needs additional 
explanation. Here one should define carefully the limiting procedure. We 
have 



1 fL/2 I +00 

^2 = J^Q2^ J ^//^ {(Ta{x)(Ta{x) + (f)a{x)(pa{x)) E ^ " (^36) 



m=—oo 
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The integral has a smooth L — oo hmit, but contains the local factor, which 
can be understood as a 5-function singularity, S{0). It is clear from eg . (11311) 
that the infinite value 

1 +CO 

- Y: 1 = ^(0) (137) 

n=— oo 

appears as a common factor on the right hand side of eq. ()13(ij) . representing 
the density of the Fourier harmonics in the interval. It must be regularized 
by cutting an upper part of the spectrum, e.g., 

1 ^ 27V + 1 

<^(0)reg = ^ E 1 = (138) 

n=—N 

where is large enough. Let us recall that the limit L — > cxd is to be taken 
afterwards. Therefore one should clarify the meaning of being large. There 
are two possibilities. One can fix 3> 1 without any relation to the size of the 
box. In this case, in the limit L oo, the (5-function vanishes. Alternatively, 
one can relate A^ with L by introducing a momentum space cutoff A: N[L) = 
LA/ {An) ^ 1. Unlike L the cutoff A has an obvious physical meaning giving 
the scale of momenta relevant for the problem. Indeed, the nth harmonic has 
a momentum pn = 27m/ L. The size of the considered box is the difference 
Pn —p-N = 4:7tN/L = A. The scale cannot be eliminated by taking the limit 
L ^ oo. One has instead 

5(0)reg = ^ + ^, (139) 

where only the second term does not contribute in the limit L ^ oo. By 
introducing the cutoff A, we suppose that the density of Fourier harmonics 
has a finite value which can be fixed phenomenologically. This is the scenario 
which we will favor in the following. 



6. The loop expansion of Z 

The perturbative series, considered in the previous section, can be resummed: 
first summing all diagrams with no closed loops (tree graphs), then those 
with one closed loop, etc. As we have already discussed, this can be done by 
generalizing the method of subsec.3.6. The tree graphs have been summed in 
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this way in ^H] , and the one-loop graphs in |2I] . The two- loop approximation 
is the next to lowest order result which contributes to the effective action, 
since the one-loop diagrams contribute solely to the amplitude. 



6.1 The two-loop approximation 

In this subsection we present a detailed computation of the two-loop approx- 
imation to the effective mesonic action generated by the functional 

q4^/-(*=l) 
-St 



Z[n, A] ~ A/'exp {i J d^xC 

+00 

J HVRa exp Q J d'xC%{n^'='^)RARB 



+00 

X 



—00 
00 

X 

n=0 



A 

E (^^ ^ABC J d'xRARBRc) (140) 



Note that in comparison with eq. ()82j) . only one critical point, related to the 
stable configuration, must be considered. It is the solution with 2 = 1 in 
eqs.dHH). We shall identify = At, and C%{n'^'=^'>) = C% in the 

following. 

By replacing the continuum of spacetime positions with a discrete lattice 
of points surrounded by separate regions of small spacetime volume Q, the 
functional integral (|140j) may be reexpressed as a Gaussian multiple integral 
over a finite number of real variables Ra{x) for a fixed spacetime point x. 
We think of VRa as the infinite product VRa —>■ Tlx ^-^^(a;), / d'^x — > Q J2x- 



+00 



Z[n, A] ~ I exp (inCst) J n di?A exp (^^nC%RARB 
^ I -00 A 

X (^1 - ^^^abc^defRaRbRcRdReRf + •••)}• (141) 

The last sum contains only terms with even powers of Ra, the odd powers 
do not contribute under the Gaussian integration. The dots mean the terms 
corresponding to the three-loop contribution and higher. 
Let us do the Gaussian integrals 

crn-r Al A/-TT I r / -o /- \ / i , ■'^'^^ ABC^ DEF^ABCDEF . . , \| 

m, A] ~ AT n I Jo exp (.^]A0 [I + ^ ,2r^^(^^2)(iV + 4) '^ + ' ' ' j j ' 

(142) 
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Here 

M = {tr C"-'y+ 6 tiC'-hiiC'-'y +8 tT[C"-')\ (143) 
and Iq is a one-loop contribution 



JV 



where are eigenvalues of the N x N matrix C^^. In our case = 18. 
C"~^ is an inverse matrix of C". The totally symmetric symbol Sabcdef 
generalizes an ordinary Kronecker delta symbol 6ab by the recurrent relation 

Sabcdef = SabScdef + Sac^bdef + . . . + SafSbcde ■ (145) 
One can find that 

27 

^ABC^DEF^ABCDEF = Q^ABC^ABC = ■ (146) 



Up to the given accuracy we have in eq. 

n ( /o e^'"''" (1 + ^F)) = n ( ^0 e^(^^-+^); 

X X 

= (n^o) (^ne*(^^-+^)j = (^n/o^ ^.^^^^=^+^70 ^^^^^ 

It shows that in the continuum limit the two-loop correction contributes 
to the effective Lagrangian as 

L-eS - L,,t + 32^^^ + 2)(iV + 4) ■ ^ ^ 

This is our final expression for the effective Lagrangian in the two-loop ap- 
proximation. 

Let us do some estimates to justify the result. For this purpose let us 
simplify the integral (jl4m) . After neglecting the symmetry group and dis- 
cretizing the spacetime it takes the form 

Z[n, A] ~ n / di?. exp {iQ (£st + ^C^^- + ^C^^) } • (149) 
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To justify the stationary phase approximation for the integral (jl49|) we as- 
sume that 

nCst > 1. (150) 

The dominating role of the Gaussian integral is reflected in the fact that 
essential values for in the integral have the order ~ 1/{QC'^^). For the 
cubic term it follows then 



where we have used that in the model considered here, ~ k, C'!.^ ~ G. If 
the parameters of the model can be chosen in such a way that the inequality 
C <^ 1 is fulfilled, the cubic power of R^ yields terms that go to zero relative 
to the Gaussian term as C — > 0, and the stationary phase approximation will 
be justified. Note that Vt may be written as an ultraviolet divergent integral 
regularized by introducing a cutoff A 

^ ' J-A/2 (27r)4 V27ry ^ ^ 

Therefore, the inequality restrics the value of A from above. 

Meanwhile, it is interesting to see that the inequality C 1 is an exact 
equivalent of ()150|) . Indeed, in the essential region, i.e., around a sharp 
minimum, one has C'^t ~ Rx^C-st, i^'it ~ Rl'^st, and so on, thus 

(O^ (Rl^st? J_ 

To summarize, the asjTiiptotical series (jl41|) with the ultraviolet cutoff 
imposed is sensible. One deals here, actually, with a series in powers of the 
dimensionless parameter (. The expansion is formally justified for ^ -C 1. 



6.2 Some final comments 

Several comments are in order: 

(1) We believe that it is the first time that loop corrections to the effec- 
tive action have been obtained at next to lowest order in the bosonization 
procedure. 
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(2) This is a nonrenormalizable theory; so we should expect an exphcit 
dependence on the cutoff in the result. This is indeed the case: the [5(0)]^ 
factor is understood as an effective arbitrary parameter (A/(27r))^ which 
together with the other coupling constants of the model, G and forms a 
dimensionless expansion parameter C, ~ k^A^/G^, the latter must be small 
to justify the loop expansion. This requirement restricts the value of A from 
above. 

(3) The field-dependent factor Ai contains all possible mesonic vertices, 
including the cr-tadpole contribution to the gap equation, and contributions 
to the masses of scalar and pseudoscalar nonets. These contributions might 
be interesting and worth studying phenomenologically. 

(4) Without future phenomenological considerations it is difficult to say 
at this stage if the loop expansion is a better approximation scheme than 
the ordinary perturbative theory considered above. It is certainly not worse 
if K is small, since the set of graphs with 2n loops includes, as a subset, all 
graphs of K^+^th order or higher in the coupling constant. Thus, in the tree 
approximation the term ~ k will dominate, the two-loop result will include 
the K, and /t^ order contributions completely, and so on. On the other hand, 
if perturbative corrections are sufficiently large, the loop expansion can be 
more appropriate. 

(5) We beheve that techniques developed here to obtain the two-loop 
contribution can be easily applied to higher orders in the loop expansion. We 
also hope that some of our findings can be used in more advanced calculations 
with non-local effective quark Lagrangians. 



7. Summary and outlook 

The purpose of this paper has been to consider several mathematical aspects 
which are related with the C/a(1) axial symmetry breaking by the 't Hooft 
determinant within the framework of the NJL model. Among them is the 
question related to the stability of the ground state, the relevance of an 
hierarchy in multi-quark interactions, and the development of techniques for 
studying the effective Lagrangian beyond lowest order in ^7a(1) breaking 
effects, so that quark - anti-quark bound-state phenomena can be examined 
in detail. 

We have shown that in this picture there is an apparent problem: the 
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model has no ground state. This conclusion is based on the stationary phase 
method applied to the generating functional of the model. This approach is a 
well established technique which allows to identify straightforwardly all criti- 
cal points, i.e., all solutions TZ^^ of the stationary phase equations associated 
with the auxiliary bosonic variables, and obtain finally the gap equations 
which are the local extrema conditions for the corresponding effective poten- 
tial. It has been shown that these gap equations have no phenomenologically 
acceptable solutions at leading order of the stationary phase approach. 

If the theory has no ground state, the theory is useless for phenomenology. 
This fact, however, can be used to obtain some constraints with a clear 
phenomenological content. Indeed, let us assume for a moment that the 't 
Hooft coupling k scales as nI . It would mean that the Lagrangian density 
£h scales with Nc, being comparable to £njl- We already know that such 
theory has no ground state and must be rejected, thus this counting for k is 
not acceptable. Indeed, all phenomenological facts related with the known 
solution of the Ua{1) problem provide for abundant evidence in favour of a 
different large Nc identification: k ~ 

There are several alternatives to save the situation. We have considered 
here the way based on the perturbative treatment of the multi-quark system. 
Two new results have been obtaind in Sec. 3 in this connection: 

(1) We have shown that the perturbative approach is related to a one 
critical point result of the stationary phase calculations: both represent the 
same function, giving different series developments for it. 

(2) We have shown that this function, if one would find it, is not equal 
to the starting integral (j2^, and is also not an asymptotical series of it. 

Both these findings are new to our knowledge and serve to understand 
better the approximations used. They probably cannot be rigorously proven 
for the functional case, but definitely solve the problem in the following 
way: by removing the destructive effect of the singular critical point in the 
generating functional, the system can be treated perturbatively around the 
stable ground state of the NJL model. Therefore, this is a reasonable step 
for qualitative phenomenological considerations at least. 

Two series expansions have been discussed in this context. We have cal- 
culated the first order corrections to both of them and gave the complete 
classification of the terms of the series, separating contributions to the effec- 
tive meson Lagrangian from the ones to the measure. We have shown that 
loop corrections serve to obtain approximately (i.e., in the framework of the 
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local model which has been considered here) the less divergent or even finite 
contributions which would originate within a more refined consideration with 
non-local multi-quark vertices. 

The new corrections, on which we report, are the direct consequence of 
the 't Hooft term in the Lagrangian and in this sense they are interesting 
by themselves and will survive even if theory and experiment disagree. Were 
this the case, it would only mean that some essential details in the multi- 
quark effective Lagrangian are still missing, thus stimulating work in this 
direction. In the opposite case, if the predictions are close to the data, our 
findings would put the model on firmer grounds. 

Another commonly used technique is the mean-field method ^Hj- It is 
easy to see that the gap equations that one gets within this approach are 
equivalent to the ones of the stationary phase method, when the latter is 
considered at leading order and the contribution of the singular critical point 
is omitted. Therefore, the effective potentials coincide in this case. It is clear 
from this comparison that the loop corrections found in our work can be also 
considered as a step beyond the leading order mean- field result. 

There is another way, which yields a rigorous mathematical solution of 
the ground state problem above: one can take into account the eight-quark or 
higher order interactions with hope to stabilize the vacuum. This approach 
will be considered elsewhere j2H]- 
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Appendix 



The Airy function 

The Airy function Ai{x) is defined as 



Ai{x) 



l/dtexp (^- + xt 



(154) 



at real values of x and can be analytically continued to the whole complex 
plane x — > z as an entire function of z. 

(a) It is real, if z is real. 

(b) It decreases exponentially for |argz| < n/S. 

(c) It increases exponentially for 7r/3 < argz < tt, and —it < aigz < — tt/S. 

(d) It oscillates at arg2; = ±7r/3,7r. 

(e) On the real axis the function has the following asymptotics 
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